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Abstract. A spectral factorization theorem is proved for polynomial rank-deficient matrix- 
functions. The theorem is used to construct paraunitary matrix-functions with first rows 
given. 
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Wiener's spectral factorization theorem [12], [4] for polynomial matrix- functions 
asserts that if 

N 

(1) S(z) = c ^ n 

n=-N 

is an m x m matrix-function (C n G <C mxm are matrix coefficients) which is positive 
definite for G T, T:= {z£ C: \z\ = 1}, then it admits a factorization 

N N 

(2) S(z) = S + (z)S-(z) = Y;AnZ n -J2 A >' n ^ z e C \W> 

n=0 n=0 

where S + is an m x m polynomial matrix-function which is nonsingular inside T, 
det S + (z) 7^ when \z\ < 1, and S~ is its adjoint, A* n = A n , n = 0,1,..., N. 
(Respectively, S~ is analytic and nonsingular outside T.) S + is unique up to a constant 
right unitary multiplier. 

The factorization (2) is also known under the name of matrix-valued Fejer-Riesz 
theorem and its simple proof is provided in [2]. Various practical applications of this 
theorem in system analysis [6] and wavelet design [1] are widely recognized. 

In the present paper we consider rank-deficient matrix polynomials and prove the 
corresponding spectral factorization theorem for them: 

Theorem 1. Let S(z) be an m x m (trigonometric) polynomial matrix-function (1) 
of order N (Cn = C*L N ^ 0) which is nonnegative definite and of rank k < m for a. a. 
z£l Then there exists a unique (up to a k x k unitary matrix right multiplier) mx k 
matrix-polynomial S + (z) = J2n=oA n z n , A n G C mxfc of order N (A N ^ 0), which is of 
full rank k for each z inside T, such that (2) holds. 

Remark. If we require of S + to be just a rational matrix-function analytic inside 
T and drop the uniqueness from the condition, then the theorem can be obtained in 
a standard algebraic manner (see [9]). Hence, as we will see below, the proof of the 
theorem provides a simple proof of the same theorem for the full rank case, k = m, as 
well. This proof is even more elementary as compared with the one given in [2] since 
it avoids an application of the Hardy space theory. 
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Prior to proving the theorem, we make some simple observations on adjoint functions 
and prove Lemma 1 on paraunitary matrix-functions. We do not claim that this lemma 
is new, but include its proof for the sake of completeness. 

If / is an analytic m x k matrix-function in C\{z 1 , z 2 , . . . , z n }, then its adjoint 

T 

f*(z) = f (1/z) is an analytic kxm matrix-function in C\{z 1; z\, ■ ■ ■ , 2*}, z* : = 
1/z, 00* = 0. Obviously, if / is analytic inside T, then /* is analytic outside T 
(including infinity). Namely, if fij(e l6 ) G Lf(T), (fa is the ijth entry of /), then 
fji(e ld ) G Lj~(T), where Lf(T) (Lj~(T)) is the set of integrable functions defined 
on T which have Fourier coefficients with negative (positive) indices equal to zero. 

T 

Since / is uniquely determined by its values on T, and f*(z) = f(z) = (f(z))* for 
\z\ = 1, usual relations for adjoint matrix-functions, like (fg)*(z) = g*(z)f*(z) and 
(f- l )*(z) = (f*)- 1 ^), etc., are valid. 

Note that if / is a rational mxm matrix-function, / G 1Z mxm , then 

(3) [f(e i6 )f*(e i6 )] .. G Loo(T) =j> fa are free of poles on T, j — 1,2,..., m, 

(Loo(T) stands for the set of bounded functions) since [f(z)f*(z)] ii = Y^j=i \fij( z )\ 2 
when \z\ = 1. 

U G TZ mxm is called paraunitary if 

(4) U(z)U*(z) = I m in the domain of U and U* , 

where I m stands for the m-dimensional unit matrix. Note that U (z) is a usual unitary 

T 

matrix for each z G T, since U*(z) = U{z) = {U{z))* when \z\ = 1 and, consequently, 

(5) U[£j T = U-\z), zeT. 

Lemma 1. If U G 7^ mxm is paraunitary and analytic inside T (its entries are free of 
poles inside T), and U^ 1 G 7^ mxm is analytic inside T as well, then U is a constant 
unitary matrix. 

Proof. The equation (4) implies that Uij(z), 1 < i,j < m, are free of poles on T 
(see (3)). Since U l3 (e id ) G L+(T) and L+(T) 9 Uj^e*) = U l3 (e i0 ) (see (5)), we have 
Uij(e %e ) G L\ (T)DL^(T). Thus Uij(z) is constant for a. a. z G T, and hence everywhere 
in the complex plane. □ 

Proof of Theorem 1. Since S is nonnegative definite on the unit circle, we have 
S*(z) = S(z), z G C\{0}. 

Observe that every polynomial matrix-function always has a constant rank in its 
domain except for a finite number of points. Without loss of generality, we can assume 
that the k x k left-upper submatrix of S, denoted by Soo, has the full rank k (a.e.) so 
that S has the block matrix form 

b[Z) -{s w (z) Su(z))> 

where S i, S w = and are matrix-functions of dimensions kx(m—k), (m—k)xk, 
and (m — k) x (m — k), respectively. Since every k + 1 rows (columns) of S(z) are 
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linearly dependent, we have 

(6) S 10 (z)S 00 1 (z)S 01 (z) = S n (z) (a.e.). 
Let 

(7) S 00 (z) = S+{z)Su(z) = S+(z)(S+)*(z) 

be the polynomial spectral factorization of S 00 which exists by virtue of the matrix- 
valued Fejer-Riesz theorem. Define 

a 10 (z) := S 10 (z)(S 00 (z)r 1 
and let So have the block matrix form 

S ° W = ©)) ' 

Then Sq(z) = [Sq (z) (Sqq(z))^ 1 S i(z)] and, taking (6) into account, one can directly 
check that 

(8) S(z) = S (z)S* (z). 

Since Sq is a polynomial matrix-function, So is a rational matrix-function, however it 
might not be analytic inside T. If is the ijth entry of Sq with a pole at a inside T, 
then we can multiply S by the unitary matrix- function U (z) = diag[l, . . . , u(z), . . . , 1], 
where u(z) — (z — a)/(l — Hz) is the jjth entry of U(z), so that the ijth entry of the 
product So(z)U(z) will not have a pole at a any longer keeping the factorization (8): 

(S U)(z)(S U)*(z) = S (z)S*(z) = S(z). 

In the same way one can remove every pole of the entries of S at points inside T. 
Thus S can be represented as a product 

(9) S(z) = S+(z)S»(z), 

where Sq is a rational matrix-function which is analytic inside T, and Sq(z) is its 
adjoint. Note that Sq(z) remains of full rank k for each z G T except possibly a finite 
number of points. 

Now, it might happen so that Sq is not of full rank k inside T everywhere. If 
\a\ < 1 and rank Sq(o) < k, then there exists a unitary matrix U such that the 
product Sq(o)U has all O's in the first column. Hence a is a zero of every en- 
try of the first column of the matrix-function Sq(z)U and the product S^(z) : = 
Sq(z)U di&g[u(z), 1, . . . , 1], where u(z) = (1 — az)/(z — a), remains analytic inside T. 
While the factorization (9) remains true replacing Sq and Sq by S* and Sf , respec- 
tively, the minors of S± will have less zeros inside T than the minors of Sq. Thus, 
continuing this process if necessary, we get the factorization 

(10) S(z) = S + (z)S-(z), 

where S + is a rational matrix-function which is analytic and of full rank k inside T. 

Now let us show that S + is in fact a polynomial matrix-function of order N. It 
suffices to show that z N ' S~{z) is analytic inside T. Indeed, since S + does not have 
poles on T (see (10) and (3)), z N S~(z) should be an analytic (on the whole C) rational 
matrix-function in this case, and therefore a polynomial. 
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It follows form (10) that 

(11) z N S-{z) = {(S+WfS+iz))- 1 ■ (S+(z)) T ■ z N S(z) 

and z N S~(z) is analytic inside T since each of the three factors on the right-hand side 
of (11) is such. 

To complete the proof of the theorem, it remains to show that the factorization (2) 
is unique, i.e. if 

S(z) = Sf(z)S^z) 

where S± is a m x k polynomial matrix-function which has the full rank k inside T, 
then 

S+( z ) = S + {z)U 

for some k x k (constant) unitary matrix U. 

Since S + (z) is of the full rank k for each z G C except for some finite number of 
singular points, there exists a matrix-function U (z) such that 

(12) S+(z) = S+{z)U{z) 
Thus U (z) can be determined by the equation 

U(z) = {(S + ) T (z)S + (z)y 1 (S + ) T (z)St(z) 

as a rational function in C. Note that U(z) is analytic inside T, and since S + and S± 
participate symmetrically in the theorem, U~ 1 {z) is analytic inside T as well. 

Due to Lemma 1, it remains to show that U G TZ kxk is a paraunitary matrix- function. 
From the equation (12), one can determine U(z) as 

U(z) = (S-(z)S + (z))- 1 S-(z)St(z) 

and, consequently, 

V(z) = St(z)S + (z)(S-(z)S + (z))- 1 . 

Hence 

U(z)U*(z) = (S-(z)S+(z))- l S-(z)Snz)-S^z)S+(z) (S-(z)S + (z)y 1 = 
= {S-(z)S + (z)y 1 S-(z)S + (z)S~(z)S+(z) (S'(z)S + (z))- 1 = I k . 

The proof of the theorem is complete. 

Remark. As one can observe, the above proof of the existence of S + is constructive. 
There are several classical algorithms to perform the factorization (7) numerically in 
the full rank case (a new efficient algorithm of such type is proposed in [5]). Further 
using the steps described in the proof, one can compute S + numerically. 

Our next theorem illustrates one of the applications of Theorem 1 in some areas of 
signal processing. Namely, mxm paraunitary matrix-functions 

N 

(13) U(Z) = Y,PnZ n = M*)] iJ= tt , Pn e C m *"\ 

n=0 

defined by (4) play an important role in the theory of wavelets and multirate filter 
banks [8] where they are known under different names, for example, lossless systems 
[11], perfect reconstruction m-filters [7], paraunitary m-channel filters [10], and so on. 



The positive integers m and iV are called the size and the length of U, respectively. 
Sometimes, the first row of a matrix-function U is called the low-pass filter, and the 
remaining rows are called the high-pass filters. Theorem 2 allows us to find the set of 
matching high-pass filters to each low-pass filter. First we give a simple proof of the 
following lemma which provides additional information about structures of paraunitary 
matrix-polynomials. 

Lemma 2. (cf. [8, Lemma 4.13]) Let (13) be a paraunitary matrix-polynomial of length 
N (p N j£ 0). Then 

(14) det U(z) = c ■ z k , where |c| = 1, andk>N. 

Proof. Since det U(z) ■ det U*(z) = 1 and det U(z) is a polynomial, it follows that 
det U (z) = cz k for some nonnegative integer k. We have 

N 1 

(15) $>>~ n = = U ~^ z ) = deU7M (C ° f = cz ' k ( CoW ^f- 

n=0 ^ ' 

Therefore k > N, since Cof U (z) is a polynomial matrix-function and p* N is not the 
zero matrix. □ 

Remark. The positive integer k in (14) is called the degree of U. Generically, a 
paraunitary matrix-polynomial U of length iV has the same degree N, although in 
some specific cases the degree is more than N. 

The following theorem was first established in [3] by a different method, however 
the presented approach gives a new insight to the problem. 

Theorem 2. For any polynomial vector-function 

(16) Ui(z) = [uu(z),u 12 (z),...,u lrn (z)], 

Uij(z) = J2n=o a jnZ n , j — 1,2, ... ,m, of length N ( Y^j=i \ a jn\ > 0) which is of unit 
norm on T 

m 

(17) \\Ui(z)\\ 2 0n = J2\u lj (z)\ 2 = l, zeT, 

there exists a unique (up to a constant left multiplier of the block matrix form ^ ^ 

where U is a (m — 1) x (m — 1) unitary matrix) paraunitary matrix-function U(z) (of 
size m and length N), with determinant cz N , \c\ = 1, whose first row is equal to (16). 

Lemma 3. Let v = (v\,v 2 , ■ ■ ■ ,v m ) T G C m be a vector of unit norm, ||v|| 2 = v*v = 
J2J=i \ v j\ 2 = 1- Then I m — vv* is a nonnegative definite matrix, 

(18) I m - VV* > 0, 

and 

(19) rank (l m — vv*) = m — 1. 
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Proof. For each column vector x G C m , we have 
x*(J m - vv*)x = ||x|| 2 



I * |2 \ II ||2 
X V > X 



*ll2n ||2 n ||2 
X V = X 



|x*|| 2 = 0. 



Hence (18) holds and x*(/ m — vv*)x = if and only if x = av for some a G C. Thus 

(19) holds as well. □ 

Proof of Theorem 2. Due to Lemma 3 and the property (17), the matrix-function 

(20) S(z) = I m -VT(z)(VT)*(z) 

is positive definite and of rank m — 1 for each z G T. (Note that the order of S is 
less than or equal to N.) Hence, by virtue of Theorem 1, there exists an m x (m — 1) 
matrix-function S + (z) of full rank m — 1, for each z inside T, such that (2) holds. 
Consequently, 



M(z) S+(z)] 



and 



U(z) 



S-(z) 

Ui(z) 
{S + ) T {z) 



is the paraunitary matrix-function we wanted to find. Indeed, clearly U(z) is of size 
m and length N, and we show that 



(21) 



detU(z) = c-z ! \ \c\ = l. 



Due to Lemma 2, det U(z) = cz k , \c\ = 1, for some positive integer k > N. Hence 
(see (15)) 



(22) 



N 

E 

n=0 



CYj n Z 



U 



z)=c-z k ■ cof (uij(z)), J = 1, 2, ... , 



m. 



Since 5* + (0) is of rank m — 1, then cof [uijifSj) ^ for at least one j G {1, 2, . . . , m} 
so that the first coefficient of the polynomial cof (uij(z)) differs from for at least one 
j. Thus it follows from (22) that k < N and hence k = N, which yields (21). The 
desired U (z) is found and let us show its uniqueness. 

Assume now that U(z) is any m x m paraunitary polynomial matrix-function, with 
the first row (16), which satisfies (21), and let U m -i(z) be the (m — 1) x m matrix- 
polynomial which is formed by deleting the first row in U (z). It is obvious that U^_ ± (z) 
is an m x (m — 1) polynomial spectral factor of (20) so that, by virtue of Theorem 

1, we get U?n-i{z) = S + (z)U <^=^ U(z) = \ \ j T ] ( f^+yr}^] immediately after we 



U \(S + ) 



establish that U m -i(z) is of full rank m — 1 for each z inside T. But rank U m -i(z) 



m — 1 for any z ^ since (21) implies that rank U(z) = m, z ^ 0, and J2n=o ® 
u\j(z) = cz~ N cof (uij(z)) (see (22)), ctjN ^ 0, implies that cof (wij(0)) ^ 0, which 
means that rank U m -i(0) = m — 1. 



jn ■ 
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